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ABSTRACT 

We study the accretion along streams from the cosmic web into high-redshift massive 
galaxies using three sets of AMR hydro-cosmological simulations. We find that the 
streams keep a roughly constant accretion rate as they penetrate into the halo centre. 
The mean accretion rate follows the mass and redshift dependence predicted for haloes 
by the EPS approximation, M oc ■ The distribution of the accretion 

rates can well be described by a sum of two Gaussians, the primary corresponding 
to “smooth inflow” and the secondary to “mergers”. The same functional form was 
already found for the distributions of specific star formation rates in observations. 
The mass fraction in the smooth component is 60 - 90%, insensitive to redshift or 
halo mass. The simulations with strong feedback show clear signs of reaccretion due 
to recycling of galactic winds. The mean accretion rate for the mergers is a factor 2 
- 3 larger than that of the smooth component. The standard deviation of the merger 
accretion rate is 0.2 - 0.3 dex, showing no trend with mass or redshift. For the smooth 
component it is 0.12 - 0.24 dex. 

Key words: methods: numerical - galaxies: evolution - galaxies: formation - galaxies: 
high redshift - intergalactic medium - cosmology: theory 


1 INTRODUCTION 

The dominant idea of galaxy form ation has changed re¬ 
cently: Decades ago it was tho ugh (iRees fc Ostrikerj|l977l : 
ISilkI Il977l : IWhite fc Reed Il978l l that galaxies collect their 
baryons through diffuse gas symmetrically falling into dark 
matter haloes and being shock-heated as it hits the gas re¬ 
siding in them - hot mode accretion. The mass of the halo 


was shown by theoretical work and simulations I 

Fardal et al. 

200ll: Birnboim & Pekell 20031: Keres et al. 

2005. 

200 c; 

Pekel & BirnboimI 

20061: lOcvirk. Pichon & Tevssierl 

2008; 

Pekel et al.l 20091. 

2 OI 3 II that galaxies at high redshift 


{z ^2), acquire their baryons primarily via cold streams of 
relatively dense and pristine gas with temperatures around 
10"* K that penetrate through the diffuse shock-heated 
medium - cold mode accretion. The peak of the stream 
activity is around redshift 3. Theoretical work predicted a 
quenching of the gas supply into h igh mass galaxies (Myjr > 
10^^ Mq) at low redshifts (z < 2) (IPekel fc Birnbohal2006ll . 
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iGenel et al.l (l201(j l showed with A-body simulations 
that about half of the dark-matter haloes’ mass is built- 
up smoothly. So as the galaxies grow bar yons should 
also be accreted semi-continuously. Indeed IPekel et al.l 
(|200ijl showed with hydrodynamical cosmological simu¬ 
lations that about two thirds of the mass are brought 
in by the rather smooth gas components, that also con¬ 
sist of mini-minor mergers with mass ratio smaller than 
1:10. The smooth and steady accretion through the cold 
streams may have been the main driver of the for¬ 
mation of the massive, clumpy and star-forming discs 
that have been observed at z ^ 2 llGenzel et al. l2008l : 
iGenel et al.ll200^ : iForster Schreiber et alT 20091 . 12011 1. Ma¬ 
jor merger events on the other hand may have only con¬ 
tributed the smaller part (|AggrtzforeYsMgLJ^AlQQ££ 1201 ll : 


ICeverino. Pekel fc Bournaiid 2010l: Geverino et al.l2015ll. 


The times scales for the gas depletion of the galaxies 
are relatively shor t compared to the Hubble time durin g all 
of cosmic history llPaddi et al.ll20lT iGenzel et al.ll201Cll l. So 
galaxies must always accrete fresh gas coming from the in¬ 
tergalactic medium in order to sustain star formation over 
such a long time at the observed level. The system is self- 
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2 Tobias Goerdt et al. 


regulated: the accretion is setting the star-f ormation rate 
indep endent of the amount of available gas dBouche et alJ 
l20inh . The galaxy’s star formation rates are fundamentally 
limited by the gas inflow. Local physics like the regula¬ 
tion by feedback seem to be only of secondary importance. 
A fact that highlights the accretion rates’ importance: the 
timescales for star formation rate and star formation both 
crucially depend on the gas accretion rate. A robust approx¬ 
imation for the ave rage growth rate of halo virial mass Myjr 
has been derived dNeistein. van den Boscli_&A3ekel| JgOOd: 

iy |l993li 


iNeistein fc Dekelll2008h using the EPS ( Lacev fc Colelll993l l 
theory of cosmological clustering into spherical haloes in 
virial equilibrium. 

Attempts to prove the cold accretion stream paradigm 
observationally are ongoing: The characteristics of Lya emis¬ 
sion produced by the cold gas streams via the release of 
gravitational energjQ have been predi cted using cosmolog - 
ical hydrodynamical AMR simulations (iGoerdt et ahll^ldl . 
The simulated Lya-blobs (LABs) can reproduce many of the 
features of the observed LABs. They can therefore be inter¬ 
preted as direct observational detections of cold stream ac¬ 
cretion. Likelihoods of observ ing streams in absor ption have 
been theoretically predicted dGoerdt et al.ll^l2h . A planar 
subgroup of satellites has been found to exist in the An¬ 
dromeda galaxy ('M31: lrbata et aO2013f l. It comprises about 
half of the population. This vast thin disk of satellites is a 
natural result of cold stream accretion dGoerdt et al.l [201^ 1 
and can therefore also be treated as observational evidence 
for the existence of cold streams. 

Observationally two main modes of star formation are 
known to control the growth of galaxies: a relatively steady 
one in disk-like galaxies and a starburst m ode which is gen¬ 
erally interpreted as driven by merging. iRodighiero et all 
d201lll quantified the relative contribution of the two modes 
in the redshift interval 1.5 < 2 < 2.5 using PACS/Herschel 
observations over the whole COSMOS and GOODS-South 
fields, in conjunction wit h previous optical/near-infrared 
data. I Sargent et al.l d2012li used their data to predict the 
shape of the infrared luminosity function at redshifts 2^2 
by introducing a double-Gaussian decomposition of the spe¬ 
cific star formation rate distribution at fixed stellar mass 
into a contribution (assumed redshift- and mass-invariant) 
from main-sequence and starburst activity. 

The characteristics of baryonic accretion on to galax¬ 
ies have been analysed with the help of cosmological hy¬ 
drodynamic simulations in various works. Key results are 
that gas accretion is mostly smooth, with mergers only be¬ 
coming important for groups and clusters. The specific rate 
of the gas accretion o n to haloes is only weak ly depen¬ 
dent on the halo mass (Ivan de Voort et~M ][2on|). The ac¬ 
creted gas is bimodal with a temperature division at 10® 
K. Cold-mode accretion dominates inflows at early times 
and declines below 2 ~ 2. Hot-mode accretion peaks near 
2 = 1 — 2 and declines gradually. Cold-mode accretion can 
fuel immediate star formation, while hot-mode accretion 

J ireferentially build s a large, hot gas reservoir in the halo 
Woods et al.ll2014h . Another third distinct accretion mode 
along with the ’cold’ and ’hot’ modes has been identified 


^ The gas remains at constant velocity as it flows down the poten - 
tial gradient towards the halo centre dGoerdt fc Ceverinoll201^ . 


dOppenheimer et al.ll201^ . called ’recycled wind mode’, in 
which the accretion comes from material that was previ¬ 
ously ejected from a galaxy. Galaxies in substantially over- 
dense environments grow predominantly by a smooth ac¬ 
cretion from cosmological filaments which dominat es the 
mass input from mergers (iRomano-Dfaz et al.ll201^ . The 
relationship between stellar mass, metallicity, and star for- 
mation rate can be explained by metal-poor gas accretion 
llAlmeida et al.l[2014l L 

The properties of cold- and hot-mode gas, are clearly 
distinguishable in the outer parts of massive haloes. The 
cold-mode gas is confined to clumpy filaments that are ap¬ 
proximately in pressure equilibrium with the hot-mode gas. 
Cold- mode gas typically has a mu ch lower metallicity and 
is mu ch more likely to be infalling (Ivan de Voort fc Schavd 
l2012li . Stream morphologies become increasingly complex at 
higher resolution, with large coherent flows revealing density 
and temperature s tructure at progressively smaller scales 
llNelson et al.l[2015h . 

The cold gas accretion rate is not an universal factor 
of the dark matter accretion rate. Galactic winds can cause 
star formation rates to deviate significantly from the exter¬ 
nal gas accretion rates, both via gas ejection and reaccre¬ 
tion. Cold accretion is broadly consistent with driving the 
bulk of the highly star-formi ng galaxies observed at 2 ~ 2 
dFaucher-Giguere et al.ll201lli . The key physical variables of 
galactic cold inflows substantially decrease with decreasing 
redshift, namely: the number of streams, the average in¬ 
flow rate per stream as well as the mean gas density in the 
streams. The stream angular momentum on t he other h and 
is found to increase with decreasing redshift llCenll20i3f . 

The average radial velocities of the inflowing material 
decrease with decreasing radius (roughly a decrease of 65% 
betw een 1 rvu and the centre of the host; I Wetzel fc Nagail 
l2015l . their figures 4 and 6, bottom panels). The velocity 
profile of the gas flowing along the streams into a galaxy’s 
halo in the form of cold streams is, contrary to what might 
be exp ected, roughly constant w ith radius instead of free- 
falling dGoerdt fc Ceveriii3l2015l j. 

In a fort hcoming companion paper dGoerdt et al.l 
lin preparatiordj we are going to address the role of merg¬ 
ers versus smooth flows by analysing the dumpiness of the 
gas streams. We will evaluate each clump mass and estimate 
a mass ratio for the expected merger. Finally we are going to 
look at the distribution of constituents amongst the clumps. 

In this paper we look at the amount of inflow - the mass 
accretion rate - as a function of radius, mass and redshift 
for the three constituents gas, stars and dark matter. The 
paper is organised as follows: in section [2] we present the 
three suites of simulations used for the analysis. In section 
Owe show individual and averaged inflow curves. In section 
Owe look at the distributions of the amount of inflow, i.e. 
at the deviations from the means presented in section O In 
section [5] we draw our conclusions. 


2 SIMULATIONS 

Galaxy snapshots of three different suites of cosmo¬ 
logical simulations based on Eulerian AMR (Adaptive 
Mesh Refinement) hydrodynamics are analysed for this 
paper. The three suites are the Adaptive Refinement 
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Streaming rates into high-z galaxies 3 



Figure 1. Simulated galaxies from two different suites of simulations (MN and ART). Shown in colour coding is the peak gas inflow along 
the line of sight. Inflow is here defined as the radial gas inflow per unit solid angle. The contours indicate gas influxes of M = 50, 
500 and 5000 Mq yr“^ sr“^, respectively. The circles refer to the virial radii. Left: one of the ART galaxies (resolution 70 pc) at 2 = 2.3, 
with Mvir = 3.5 X IO^^Mq. Right: a typical MN galaxy (resolution Ikpc) at 2 = 2.5, with Mvir = IO^^Mq. The inflow is dominated in 
both cases by three cold narrow streams that are partly clumpy. 


Tree (hereafter ART; 

Ceverino. Dekel & Bournaud 20101: 

Ceverino et al.l I2OI2I: 

Dekel et al. 2013: Ceverino et al.l 

20151). the Arp ICeverino et al.ll2014l; IZolotov et al.ll2015l) 


as well as the Horizon-Mare Nostrum (hereafter MN; 
lOcvirk. Pichon fc Tevssi^l2008l ) simulation. The ART and 
the Arp suites of simulations both consist of several simu¬ 
lations zooming in with a maximum resolution of 15 — 70 pc 
at 2 = 2 on individual galaxies that reside in dark-matter 
haloes of masses (0.07 — 1.98) x IO^^Mq at 2 = 2.3. The 
MN simulation contains hundreds of massive galaxies in a 
cosmological box of side 71 Mpc with a maximum resolution 
of ~ 1 kpc. 

To show the morphology of the gas streams we present 
in figure [T] gas influx maps of two sample galaxies from the 
MN and the ART suite. The panels demonstrate the dom¬ 
inance of typically three, narrow cold streams, which come 
from well outside the virial radius along the dark-matter fil¬ 
aments of the cosmic web, and penetrate into the discs at 
the halo centres. The streams are partly clumpy and partly 
smooth, even in the simulations with higher resolution. 


2.1 High-resolution ART simulations 

The ART simulations were run with the AMR code ART 
(Adaptive Refinemen t Tree; iKravtsov. Klvpin fc KhokhlovI 
I1997I : lKravtso^l2003l ) with a spatial resolution better than 
70 pc in physical units. It incorporates the relevant phys¬ 
ical processes for galaxy formation that are: gas cooling, 
photoionisation heatin g, star formation, metal e nrichment 
and stellar feedback JCeverino &: Klvoinl l2009l) . Cooling 
rates were computed for the given gas density, tempera- 
ture, metallici t y, an d UV background based on CLOUDY 
llFerland et al.l Il998l ) . Cooling is assumed at the centre 
of a cloud of thickness 1 k pc llCeverino- Rodrigu^ I2OO8I : 
iRobertson fc Kravtsovl[2008l ) . Metallicity dependent, metal¬ 
line cooling is included, assuming a relative abundance of 
elements equal to the solar composition. The code imple¬ 
ments a “constant” feedback model, in which the combined 


energy from stellar winds and supernova explosions is re¬ 
leased as a constant heating rate over 40 Myr (the typical 
age of the lightest star that can still explode in a type-II 
supernova). Photo-heating is also taken into account self- 
consistently with ra diative cooling. A unifor m UV back¬ 
ground based on the lHaardt fc MadarJ (|l99fil ) model is as¬ 
sumed. Local sources are ignored. In order to mimic the self¬ 
shielding of dense, galactic neutral hydrogen from the cos¬ 
mological UV background, the simulation assumes for the 
gas at total densities above n = 0.1 cm“® a substantially 
suppressed UV background (5.9 x 10^® erg s“^ cm“^ Hz“^, 
the value of the pre-reionisation UV background at 2 = 8). 
The self-shielding threshold of n = 0.1 cm“^ is set by ra¬ 
diative transfer simulations of ionising photons coming from 
the cosmological bac kground as well as from the galaxy itself 
(see appendix A2 of iFumagalli et ^l201ll ). 


The ART code has a unique feature for the purpose of 
simulating the detailed structure of the streams. It allows 
gas cooling to well below lO^K. This enables high densi¬ 
ties in pressure equilibrium with the hotter and more di¬ 
lute medium. A non-thermal pressure floor has been imple¬ 
mented to ensure that the Jeans length is resolved by at least 
seven resolution elements and thus p revent artificial frag¬ 
mentation on the smallest g ri d scale llTruelove et al. 1997; 


R.obertson fc Kravtsov! I2OOS ; ICeverino. Dekel fc BournaudI 


20inl) . It is effective in the dense (n > 10 cm ®) and cold 


(T < ImK) regions inside galactic disks. 


The equation of state remains unchanged at all den¬ 
sities. Stars form in cells where the gas temperature is 
below lO'^K and the gas density is above a threshold of 
n = lcm“® ac cording to a stoc hastic model that is con¬ 
sistent with the iKennicrittl (Il998l ) law. In the creation of a 
single stellar particle, a constant fraction (1/3) of the gas 
mass wit hin a cell is converted in to stellar mass (see ap¬ 
pendix of ICeverino fc KlvDinll200^ ). The ISM is enriched by 
metals from supernovae type II and type la. Metals are re¬ 
leased from each star particl e by SNII at a consta nt rate 
for 40 Myr after its birth. A iMiller fc Scald (Il979l ) initial 
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4 Tobias Goerdt et al. 


mass function is as s umed which is matching the results of 
IWooslev fc Weaver! (Il995l l. The metal ejection by SNIa as¬ 
sumes an exponentially declining SNIa rate from a max¬ 
imum at 1 Gyr. The code treats the advection of metals 
self-co nsistently and it distinguish es between SNII and SNIa 
ejecta lICeverino-R odrigue jl2nn8f l . 

The initial conditions for the ART simulations were 
created using low-resolution cosmological A^-body simula¬ 
tions in comoving boxes of side 29 - 114 Mp c. Its cosmologi¬ 
cal p arameters were motivated by WMAP5 (iKomatsu et alJ 
1200^ 1. The values are: = 0.27, Ha = 0.73, fib = 0.045, 

h = 0.7 and erg = 0.82. We selected 31 haloes of Mvir — 
10^^ Mq But z = 1.0. For each halo, a concentric sphere of ra¬ 
dius twice the virial radius was identified for re-simulation 
with high resolution. Gas was added to the box following 
the dark matter distribution with a fraction /b = 0.15. The 
whole box was then re-simulated, with refined resolution in 
the selected volume about the respective galaxy. The dark 
matter particle mass is 5.5 x lO^M©, the minimum star par¬ 
ticle mass is 1O^M0, the smallest cell size is 35 pc (physical 
units) at z = 2. 

2.1.1 Simulations including radiation pressure (Arcp) 

The Arp suite of simulations llGeverino et al.l l2014l : 
IZolotov et al.l 120151 1 is a further development of last sub¬ 
section’s ART suite: Apart from the features already pre¬ 
sented there, it also includes the effects of radiation pres¬ 
sure by massive stars. The radiation pressure was modelled 
as a non-thermal pressure that acts only in dense and op¬ 
tically thick star-forming regions in a way that the ionis¬ 
ing radiation injects momentu m around massive s tars, pres¬ 
surising star-forming regions l|Agertz et ahllioi^ . their ap¬ 
pendix B). The Arp suite uses a stronger feedback that 
brings the stellar-to-halo mass ratio closer to the estimates 
by ab unda nce matching (iMoste r et al.lj2 013l: iBehroozi e t al.l 
l2013l ~). See lGeverino et ahl ll2014l l and iMoodv et an ~i 2014h . 
The adaptive comoving mesh has been refined in the dense 
regions to cells of minimum size between 17 — 35 pc in phys¬ 
ical units. The DM particle mass is 8.3 x 10^ M©. The parti¬ 
cles representing star clusters have a minimum mass of 10® 
Mq , similar to the stellar mass of an Orion-like star cluster. 
The initial conditions for the Arp simulations were created 
using low-resolution cosmological A-body simulations in co¬ 
moving boxes of side 14 - 57 Mpc. We selected 34 haloes of 
Mvir = 1.18 - 14.7 X 10“Mo at z = 1.0. 


2.2 RAMSES Horizon-MareNostrum simulation 

The MN simulation uses the AMR code RAMSES llTevssieil 
l2002l 'l. The spatial resolution is ~ 1 kpc in physical units. 
UV h eating is included assuming the iHaardt fc Madaul 
lll996lj background model, as in the ART simulation. The 
code incorporates a simple model of supernovae feed¬ 
back and metal enrichment u si ng th e implementation de¬ 
scribed in iDubois fc Tevssi^ (120081 '). The cooling rates 
are calculated assuming ionisation equilibrium for H and 
He, inc luding both collisional ionisation and photoioni¬ 
sation (iKatz. Hernguist fc Weinberd [l992 l. Metal cooling 
is also included using tabulated CLOUDY rates, and is 
assumed proportional to the metallicity, relative to the 


label 

suite 

Mvir [IO12 Mq] 

Z 

Agal 

IQl® 

MN 

10.47 ± 0.56 

1.57 

12 

IQl® 

MN 

10.49 ± 0.93 

2.46 

12 

5 X I012 

MN 

5.00 ± 0.045 

1.57 

12 

5 X IOI2 

MN 

5.48 ± 0.26 

2.46 

11 

I012 

MN 

1.03 ± 0.003 

1.57 

8 

IOI2 

MN 

1.01 ± 0.004 

2.46 

12 

I012 

MN 

1.03 ± 0.006 

4.01 

9 

IQll 

MN 

0.099 ± 0.000 

1.57 

12 

10“ 

MN 

0.099 ± 0.000 

2.46 

7 

10“ 

MN 

0.099 ± 0.000 

4.01 

12 

1.9 X 10“ 

ART 

1.907 ± 0.217 

1.14 ± 0.02 

34 

1.3 X 10“ 

ART 

1.286 ± 0.093 

1.60 ± 0.02 

73 

8.6 X 10“ 

ART 

0.863 ± 0.046 

2.25 ± 0.02 

109 

3.9 X 10“ 

ART 

0.391 ± 0.034 

3.40 ± 0.04 

119 

7.1 X 10“ 

Arp 

0.707 ± 0.055 

1.14 ± 0.02 

41 

6.7 X 10“ 

Arp 

0.672 ± 0.049 

1.58 ± 0.02 

47 

4.9 X 10“ 

Arp 

0.491 ± 0.028 

2.21 ± 0.03 

57 

2.9 X 10“ 

Arp 

0.290 ± 0.019 

3.27 ± 0.05 

62 

2.6 X 10“ 

Arp 

0.260 ± 0.010 

1.13 ± 0.02 

40 

2.4 X 10“ 

Arp 

0.242 ± 0.010 

1.59 ± 0.02 

50 

1.6 X 10“ 

Arp 

0.160 ± 0.007 

2.29 ± 0.03 

60 

7.3 X IQi® 

Arp 

0.073 ± 0.004 

3.53 ± 0.05 

66 


Table 1. The various bins of galaxies for our analyses. They are 
used throughout the whole paper. ’Label’ is the tag the corre¬ 
sponding bin is labelled by in our figures. It is a mass close to the 
ensemble’s actual mean virial mass. ’Suite’ denotes the suite of 
simulations the bin stems from. Mvir is the mean virial mass of the 
bin together with its standard deviation, z is the mean redshift of 
the ensemble together with its standard deviation. Galaxies from 
different MN snapshots are not combined therefore the redshift’s 
standard deviation is always zero for all MN bins. Ag^i gives the 
total number of galaxies in the given bin. 

iGrevesse fc Sau\^ lll998ll solar abundances. Unlike in the 
ART simulation, no cooling below T < 10‘^K is computed, 
and no self-shielding of the UV flux is assumed. 

For high-density regions, the RAMSES code considers a 
polytropic equation of state with 70 = 5/3 to model the com¬ 
plex, multi-phase and turbulent struc t ure of the inter-stellar 
medi um (ISM) (lYeoes et al. 1 19971: ISpringel k, HernauistJ 
200j) in a simplifie d form (see Schave fc Dalla Vecchiall2008l : 
Dubois fc Tevssi^ l2008l l. The ISM is defined as gas with 
hydrogen density greater than tir = 0.1 cm“®, one order of 
magnitude lower than in the ART simulation. Star forma¬ 
tion has been included, for ISM gas only, by spawning star 
particles at a rate consistent with the iKennicuta (Il998l l law 
derived from local observations of star forming galaxies. 

The MN simulation implemented a pressure floor in 
order to prevent artificial fragmentation, by keeping the 
Jeans lengthscale, Aj oc larger than the size of 

four grid cells everywhere. At every position where n > 
0.1 cm“®, a density dependent temperature floor was im¬ 
posed. It mimics the average thermal and turbulent pres¬ 
sure of the multiphase ISM dSoringel fc Hernauisd l2003l : 
IPalla Vecchia fc Schavd I20n8h . The gas is allowed to heat 
up above this temperature floor and cool back again. The 
temperature floor follows a polytropic equation of state with 
Tfloor = To(n/no)®'°~®, where To = 10"* K and no = 0.1 
cm~®. The resulting pressure floor is given by Pfloor = 

n-H fee Tfl oor • 
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For each stellar population, 10% of the mass is assumed 
to turn into supernovae type II after 10 Myr, where the 
energy and metals are released in a single impulse. For each 
supernova, 10% of the ejected mass is assumed to be pure 
metals, with the remaining 90% keeping the metallicity of 
the star at birth. SNIa feedback has not been considered. 
The initial conditions of the MN simulation were constructed 
assuming a ACDM universe with Hm = 0.3, Ha = 0.7, Hb = 
0.045, h = 0.7 and as = 0.9 in a periodic box of 71 Mpc. 
The adaptive-resolution rules in this simulation were the 
same everywhere, with no zoom-in resimulation of individual 
galaxies. The dark matter particle mass is 1.16 x IO^Mq, the 
star particle mass is 2.05 x 10 ®Mq, the smallest cell size is 
1.09 kpc physical, and the force softening length is 1.65 kpc. 

We will always show averaged results for an ensemble of 
galaxies having very similar masses and redshifts. In table 
[T] we show a summary of the various bins of galaxies we 
use. Since we have better statistics for the MN simulation 
we were able to bin galaxies from a narrower mass range 
therefore the standard deviations of the mean mass of a MN 
bin is usually much smaller. To partly compensate for that 
we combine ART galaxies from adjacent redshifts increasing 
the statistics but introducing a standard deviation into the 
mean redshift of the bin. 


3 AMOUNT OF INFLOW 

First the average inflow is computed as a function of ra¬ 
dius from the simulations. To do so the amount of mass 
is measured that crosses within a small time At a spheri¬ 
cal shell of radius r that is centred around any given host 
galaxy. This is done independently for gas, stars and dark 
matter, which are the three different constituents in the sim¬ 
ulations. To get an inflow rate, we divide the mass that is 
crossing through a given shell by the time At taken. The 
reader should note that we do include the expansion due to 
the Hubble flow in this calculation. As with most cosmolog¬ 
ical codes only the peculiar velocities are written out into 
the snapshots onto hard disk. We then add the velocity of 
the Hubble expansion as measured from the centre of the re¬ 
spective host halo by hand. The crucial difference between 
stars and dark matter particles on the one hand and gas on 
the other hand is that stars and dark matter are collisionless 
and therefore virialised within the virial radius Cvir, whereas 
gas on the other hand is neither. In case of the gas inflow 
we chose to only take into account inwards radial velocity 
cells in order to get the cleanest estimate for the amount of 
the inflow. Whereas for stars and dark matter we have to 
take into account all the material (i.e. only the net influx). 
If one takes into account only the inflowing stellar particles, 
the difference would not be very big, the stellar inflow rates 
would only be ~ 25% bigger. The stars follow the gas very 
precisely. For the dark matter particles on the other hand 
the differences are tremendous: The inflow rates increase by 
more than a factor of two. It is worth stressing that we do 
not use any information about the temperature: so hot in¬ 
flowing gas is taken into account, too. The hot gas however 
contributes very little to the inflow, since the vast majority 
of the mass of the material is in the cold phase. 

We also want to briefly look at the inflow rates of all 
cells to compare those to the ones from our choice of taking 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

r [Cvir] 


Figure 2. The average inflow as a function of radius together with 
its 1—(T standard deviations. The MN (dotted blue), the ART 
(solid red) and the Aj^p (dashed green) simulations are compared. 
In gas (top panel) an average “ground” inflow is seen in gas which 
is fairly constant over the different radii. On top of this smooth 
ground inflow there are distinct peaks in the inflow, corresponding 
to bigger clumps. In the stellar distribution (middle panel) signals 
of inflow are only seen at those gas peaks. Stars seem to form 
only within those gas clumps. The stellar inflow rates are noisy 
but roughly constant. In dark matter (bottom panel) the inflow 
increases with radius. The gas inflow of the App simulation is 
fairly high due to its recycling processes. 
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into account only the cells which have an inwards radial ve¬ 
locity. The results vary for the different suites of simulations: 
for the MN suite, the one with the weakest feedback, there 
is hardly any difference: the inflow rates are still constant 
at roughly the same level. For the ART simulations which 
has stronger feedback implemented, the average inflow rate 
drops by a factor of two, when also taking into account cells 
with outflowing velocities. Also we see that the inflow rates 
are smaller toward s the centre of the host halo. A n effect al¬ 
ready described by Ivan de Voort fc Scha*^ (I 2 OI 2 I . when tak¬ 
ing into account only inflowing gas, they find approximately 
constant accretion rates as well). For the Arp suite of simu¬ 
lations which has by far the strongest feedback the average 
inflow rate drops by a factor of ten, but for this suite the in¬ 
flow rates are hi gher towards the centre of the h ost halo. We 
see an inverted Ivan de Voort fc Schavd (I 2 OI 2 I I effect. The 
differences between taking outflowing gas cells into account 
or not might seem large, the reader however should note 
that we are mainly interested in the inflowing material and 
not so much the outflowing material, so concentrating on 
the inflowing cells only is the right choice. 

We compute inflow rates for the galaxies. For reduc¬ 
ing the statistical noise that is present in the inflow of a 
single galaxy we stack the amount of inflow of all available 
galaxies having similar redshifts and masses, from all three 
suites of simulations galaxies that are available at a variety 
of different halo masses and redshifts. Since the MN is a 
fully cosmological simulation it is possible to compile bins 
of galaxies spanning several orders of magnitude in mass for 
various redshifts. The ART suite consists only of resimulated 
galaxies at comparable masses and therefore there is at any 
given redshift only one mass bin available which evolves in 
mass with cosmic time. The situation of the Arp suite is 
similar, however for this suite two different mass bins are 
available at the given redshifts, but those masses also evolve 
with cosmic time. The average halo masses and redshifts of 
the chosen bins are summarised in table [T] 

In figure [ 2 ] the average inflow as a function of radius for 
galaxies at z = 2.2 — 2.5 with Mvr = 5 — 10 x 10^^ Mq of all 
three MN, ART and Arp are shown. The galaxies have an 
average “ground” inflow in gas of about 100 Mq yr“^ which 
is fairly constant over the different radii. This is the smooth 
component of the inflow, the cold streams. On top of this 
smooth ground inflow there are a couple of distinct peaks. 
These peaks can still be seen in the averaged data, especially 
if the number of galaxies in a certain bin is small (i.e. the 
MN simulation). Those peaks correspond to bigger clumps 
which are the merger events. We use the term “merger” to 
describe any major or minor merger of mass ratio pm -5 0.1, 
as distinct from “smooth flows”, which include “mini-minor” 
mergers with mass ratios /im < 0.1. We will refer to this 
bimodality of the modes of inflow, smooth accretion versus 
clumps or merger events, throughout the whole paper. 

The MN simulation is in general more noisy which can 
be attributed to its lower resolution. The Arp simulation 
has a lower stellar and dark matter infall because of its much 
lower virial mass. The gas infall of the Arp simulation on 
the other hand is higher than that of the ART or MN simu¬ 
lation, which is due to the fact that the Arp simulation has 
stronger feedback than the ART suite by design. Addition¬ 
ally to all the feedback, the ART suite of simulation already 
has, the Arp simulation also includes feedback by radiation 


pressure. ICeverino et al. liH) compared the effects of the 
feedback models used in ART and Arp suites of simula¬ 
tions. SFR and stellar mas ses were lower by a factor ~ 3 
with the stronger feedback. IZolotov et al.l (I2OI5I I computed 
the outflow rates (their figures 2 and 3) of the Arp suite 
and they found outflow rates similar or even higher than 
the SFR. Therefore the Arp suite has more massive out¬ 
flows and part of the outflowing material can rain down back 
to the central galaxy. An effect usually coined “recycling” 
JOppenheimer et al.(l201^ ') . As we will see later in figure [3] 
the gas inflow in the Arp suite of simulations is enhanced 
due to recycling compared to the MN or the ART simula¬ 
tions. The inflow of the Arp simulations is higher by a factor 
/rec(Af) = 4.0 compared to the MN or the ART simulations 
due to this effect. This factor for the inflow stays constant 
over the whole redshift range and also over the whole host 
halo mass range considered. We will also see a similar effect 
of recycling later in the paper during the fitting procedure 
of the distributions. 

Otherwise all three suites behave similar: The gas inflow 
is roughly constant at around 100 Mq yr“ MTeklul (l2012li 
reports a similar behaviour in SPH simulations. She finds a 
total gas inflow M = 90 Mq yr“^ for a galaxy with Mvir = 
1.5 X 10*^^ Mq at z = 2.33. Coming back to our own analysis 
and figure [ 2 I we see that the MN simulation seems to have a 
marginally higher averaged gas inflow than the ART which 
can be accounted for by the 16% more mass the average 
MN galaxy has. The gas inflow of the Arp simulation in 
turn is even higher, which can be accounted for by either 
recycling or more efficient cooling from more efficient metal 
enrichment. 

In MN we see that there is less stellar infall, but more 
gas infall by roughly a factor of ~ 2 compared to the ART 
simulations. In Arp there is even less stellar infall and even 
more gas infall. The stellar inflow behaviour in general is 
extremely noisy. Its inflow rates oscillate between a few and 
up to about 80 Mq yr“^ for both, with the ART and the 
Arp being more constant compared to MN, which is due to 
the fact that the bins of the ART and the Arp simulations 
consists of considerably more galaxies than the bins of the 
MN (see table [T}. Any peaks in the inflow profile (e.g. due 
to merger events) are much more smoothed out by simply 
the averaging process. It is noteworthy that in the stellar 
distribution there is only signal of inflow at the positions of 
the very peaks of the gas inflow, indicating that the smooth 
stream inflow component is fairly depleted from stars. The 
stars seem to flow into the central galaxy in the form of 
bigger clumps only. It is in this sense that the stars seem 
to follow the positions of the gas as we will refer to later. 
The lines for the dark matter inflow behave similar at larger 
radii. They exhibit slightly decreasing infall between 400 Mq 
yr“^ at 1.7 r^ir down to 200 Mq yr“^ at 0.7 Cvir- Inside of 
~ 0.7 rvir the three suites differ: MN stays constant although 
the scatter is larger, whereas ART and Arp continue to 
decrease. The ART decreases down to 50 Mq yr“^ at 0.2 
Tvir and the Arp down to 50 Mq yr“^ already at 0.8 rvir- 

The averaged amount of inflow into the galaxies as a 
function of radius is by and large constant with radius r 
in all of the mass-redshift bins. This does not come as a 
surprise since the inflowing mass should be roughly con¬ 
served. The following systematic variations of the amount 
of inflow as a function of host halo mass Mvir and redshift 
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Figure 3. The amount of inflow of the gas as a function of halo mass M^ir (left panel) or as a function of redshift (right panel). 

The ART data is plotted as filled squares, the MN data as open circles and the Aj^p data as filled triangles. The data points have 
been rescaled according to (M) oc [1 + 2 ]^ ® (left panel) or (M) oc (right panel) in order to match 2 = 2.46 (left panel) or 

A^vir = 10^^ Mq (right panel). The original redshift or host halo mass values are given in colour coding. The App values are higher 
than the other two simulations, but they also follow clearly the same slopes (indicated by the dashed grey line). This is due to the fact 
that the App simulations have more re-accretion due to recycling. One should corre ct the gas inflow i n tho se simulations by a factor 
of fiec{M) = 4.0 due to this effect. The solid black lines are the predictions from the iNeistein fc Dekell <l2008h model (equation [2} . The 
agreement between the data and the model is truly remarkable (Taking /rec(M) into account for App). 


2 become apparent: (a) The averaged amount of inflow in¬ 
creases with increasing host halo mass, (b) Both, the amount 
of inflow as well as the halo mass increase by roughly the 
same factor (one order of magnitude higher halo mass re¬ 
sults in roughly one order of magnitude higher inflow) and 
(c) the total amount of inflow decreases slightly with in¬ 
creasing time (decreasing redshift). These trends are most 
convincingly seen in gas. 

In the following we will make frequent use of basic sta¬ 
tistical measures of our samples, namely the linearly aver¬ 
aged mean (M), the mean of the logarithm (/tq) and the 
standard deviation (cro). The reader should note that those 
(unlike the later defined /ri, fj, 2 , ni, (J 2 or ai/2) are intrinsic 
properties of the distributions, that are independent from 
any model, equation or fits to it. Their definitions are: 


1 

(M) = -J2m, 

i=l 
1 """ 


MO 


CTO = 


1 

\ " i=l 


( 1 ) 


Us is the total number of shells in the sample, so the number 
of galaxies in the respective bin times the number of different 


radii which are used. Mi is the amount of gas inflow through 

the respective shell. _ 

The theoretical prediction llNeistein fc Dek^l2008h of 
the average gas inflow rate (M) into a host halo as a function 
of virial mass Mvir and redshift 2 is given by 


(M) = 100 — 
yr 



( 2 ) 


To compare our results to this prediction, the magnitude 
of the inflow in the simulations will be averaged over the 
whole radius range for all mass-redshift bins. The result¬ 
ing averaged inflow values Minflow of the gas are shown in 
figure |3] as a function of halo mass Mvir and redshift 2. In 
order to match 2 = 2.46 or Mvir = 10^^ M0 in figure [3] the 
data points had to be re-scaled according to the above equa¬ 
tion (i.e. the data in the left panel was rescaled according to 
(M) oc [l-|-2]^ ® and the data in the right panel was rescaled 
according to (M) oc [Mvir]^'^®)- The original values of red¬ 
shift (left panel) or host halo mass (right panel) are given in 
colour coding. Equation © is overlain as solid black lines 
in both panels. 

The simulated values for the inflow of the MN and the 
ART simulations are consistent with the theoretical predic¬ 
tions. The agreement is robust over three orders of magni¬ 
tude in halo mass and also over the entire redshift range. 
The simulated inflow values of Arp on the other hand are 
higher than the other two suites of simulations and the the¬ 
oretical prediction. Interestingly the Arp suite clearly fol- 
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logio(A^ [Mq yr l]) logio(j'^ [M© yr l]) logio(M [M© yr l]) 


Figure 4. Normalised distributions of the gas inflow of selected mass-redshift bins for MN (left panel), ART (centre panel) and Arp 
(right panel) in the radius range from 0.4 — 1.7 r^vir- The solid blue histograms are the real distributions from the simulations. Shown 
in red are the individual double-Gaussian (equation [3jl fits to that particular panel. A double-Gaussian consists of the sum of a first 
Gaussian (plotted in grey) and a second Gaussian (plotted in green), po and (jq denote the statistical measures of the real underlying 
distribution, whereas pi, p 2 , cri, cr 2 and ai /2 from equation @ are the best fit parameters for the double-Gaussian. In all of these panels 
there is again a very beautiful double-Gaussian behaviour, the quality of the fits is quite striking in all three panels. It is most impressive 
for the right panel, since the Arp suite is our technically most advanced simulation. 


lows the same slopes with respect to both, the redshift and 
the halo mass, still. This parallel shift by a constant factor 
/rec(M) = 4.0 is due to the fact that the Arp simulations 
have more massive outflows and therefore more reaccretion 
due to recycling. Since equation m describes only the pri¬ 
mary accretion via cold streams and not the secondary ac¬ 
cretion due to recycling, the Arp simulations should be 
scaled down by this factor. We will also see a similar effect 
of recycling later in the paper during the fitting procedure 
of the distributions. 

Accounting for recycling in the way described leads to 
an excellent agreement between the theoretical prediction 
and all the simulations including the Arp suite, too. Only 
the highest redshift / lowest mass bins of both the MN and 
the ART suite of simulation lie a little bit further away from 
the relation. This is most probably due to the fact that these 
two bins are in a point in parameter space where the simu¬ 
lations have a lower resolution and might therefore be more 
vulnerable to numerical effects. 

Gas supply of high m ass galaxies (Myir > 10 ^^ M©) 
is theoretically predicted dPekel fc BirnboimI l2006tl to be 
quenched at low redshifts (^ < 2). It would be nice to spot 
this effect in simulations. However this turns out to be a 
difficult task to do: The ART galaxies are all at roughly 
Mvir ~ 10^^ M©, the Arp galaxies are even lighter and 
both are therefore too light to test this prediction. Neither 
the 10^^ M© nor the 5 x 10^^ M© nor the 10^^ M© bins 
of the MN simulations (table [T]) show drops of the gas in¬ 
fall from 2 = 2.46 down to z = 1.57 (see the values pre¬ 
sented in figure (Sj that exceed the statistical uncertainties 
(compare the 1-cr standard deviations shown in figure [2]). 
Unfortunately the simulations do not produce galaxies in 
such high mass regimes at any earlier epochs sufficiently 
frequent enough. More high resolution simulations of those 
high sigma peak galaxies are needed to establish a possi¬ 
ble quenching of gas supply at low redshifts for high mass 


haloes in simulations. However such simulations would be 
computationally extremely demanding. 


4 DISTRIBUTIONS 

In section O the average amount of inflow was presented. In 
this section, the behaviour of the underlying distributions is 
discussed as a function of host halo mass and redshift. So the 
likelihoods of certain deviations from the means presented 
in section [3] are determined. For this the amount of inflow 
is measured through thin spherical shells at a number (200) 
of different radii between 0.4 — 1.7 Vvii around all available 
galaxies separately. This is done for every redshift and mass 
bin from all three suites of simulations. Having collected 
all these data no average as in section [3] is computed, but 
instead histograms how often a certain inflow value occurs 
are compiled. All measurements of the amount of inflow are 
treated equal, regardless of the radius they were obtained at. 
The resulting normalised distributions of the gas inflow are 
plotted in figure [4] for all three suites of simulations (MN, 
ART as well as Arp) as solid blue histograms. Similar figures 
showing the distributions of the baryonic inflow instead of 
the gas were produced, found to look extremely similar and 
therefore omitted. Some trends are seen, most striking: the 
overall mean {no, as defined in equation [T]) of the sample 
distributions (examples are shown in figure [4] as solid blue 
histograms) appears to increase with increasing mass and 
also with increasing redshift. 

The best description for those histograms is the sum 
of a strong first Gaussian (plotted in grey) with a smaller 
mean {fii) and an also smaller standard deviation (ai) and 
a weaker second Gaussian (plotted in green) having a bigger 
mean {^ 2 ) and a bigger standard deviation (0-2). This sum 
(plotted in red) can be expressed as the following “double- 
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Figure 5. Results of the fits from equation ((S]! to the gas distributions of all of the bins mentioned in tabled examples are plotted in 
figure [4| Shown are the variations of the five parameters /ii, ^ 2 , crit (^2 and ai /2 from equation (O as a function of halo mass Mvir and 
redshift z. These are the parameters of a fit. In this plot the data is completely unrescaled, unlike the same data shown in figure 0 All 
Arp values for /ii and ^2 are due to effects of recycling higher than expected. Physically these parameters show that (a) there is more 
absolute inflow into heavier haloes, (b) there is more absolute inflow at earlier cosmic times and there are heavier mergers happening 
(c) at earlier epochs and also for (d) higher mass host haloes. An analogue plot for the baryons was also produced, since it looks very 
similar to this one it was omitted. 


Gaussian” distribution: 

-{x- 

2a\ 

1 — ai/2 ~{x — r-2)^ 

CT2V^ [ 2(T2 


fix) = 


0 - 1/2 


01 


exp 


(3) 


The variable ai/2 gives the fractional strength of the first 
to the second Gaussian. This model reflects the two dif¬ 
ferent physical origins of the inflow: The majority of the 
gas inflow is expected to flow into the host halo by smooth 
accretion and a much smaller amount is expected to en¬ 
ter the host by merger events. Each of these mechanisms is 
accounted for by one of the two Gaussians of the distribu¬ 
tion. Inflow by smooth accretion will naturally occur with 
a lower amount of inflow at any one instant (small fj,i) but 
contribute much more to the inflow (big ai/2) whereas inflow 
by merger events will have a much higher amount of inflow 
at any one instant (big fi 2 ) but contributes much less to the 
inflow (small 1 —ai/2). So the first Gaussian of this distribu¬ 
tion represents the infall by smooth accretion whereas the 
second Gaussian of this distribution represents the infall by 
merger events. 


It is useful to introduce the double Gaussian functional 
form at this point for several reasons: first, because it is 
a neat and very accurate description of the behaviour of 
the inflow distribution in hydrodynamical simulations. Sec¬ 
ondly it reflects nicely the physics of the bimodal nature 
of the inflowing material (mergers vs. smooth inflow). Last 
but not least this functional form’s parameters teach us di¬ 
rectly the mean accretion rate for the mergers (fi 2 ) or for 
the smooth component (/ii), their standard deviations (a 2 , 
(Ti) and maybe most importantly the relative strength of 
accretion via the smooth component compared to the accre¬ 
tion via merger events (01/2). The above formula has already 
been used to describe the shape of the specific star formation 
rate distribution at redshifts 2: ^ 2 in observational work 
llSargent et al.l[2012l . their equation 1). The similarity indi¬ 
cates that smooth accretion might fuel main-sequence star 
formation on the one hand and accretion through merger 
events might fuel star formation in starburst events on the 
other hand. 

The overall equation (shown in red) is an excellent de¬ 
scription for all three suites of simulations. Only exceptions 
are the very low mass (Mvir = 10^^ M©) MN galaxies which 
are noisy and the low redshift (z ^ 1.7) ART cases in which 


© 2015 RAS, MNRAS OOO.fTHTSl 
















10 Tobias Goerdt et al. 


Z + l ^vir [M©] 

2.5 3 3.5 4 4.5 5 10“ 


IQii I0i3 2.5 3 3.5 4 4.5 5 



Figure 6. Shown are the positions of the centres of the first upper panels) and the second {fi 2 , lower panels) Gaussians of the 
distributions of the inflow as a function of mass (left panels) or as a function of redshift (right panels). There are very strong logarithmic 
dependencies present in all panels as described by equation m which is overlain as solid black lines. The data points within the left 
panels were scaled with the help of the scaling relation in the corresponding right panel to the values indicated. The colour bar axes 
indicate the values the data points used to have before rescaling. Please refer to figure [5] for the unrescaled data. The panels show that 
our simple phenomenological model (equation and table [2]| is an excellent description of the data. 


either the second Gaussian (green) is tiny or its mean ( 112 ) 
is not much bigger than the mean (^ 1 ) of the first Gaussian 
(grey). The double-Gaussian behaviour can most beautifully 
been seen in the technically most advanced Arp suite of 
simulations. Here all bins show an almost prototypical be¬ 
haviour in that sense. 

We now want to look at the characteristics of the model. 
Therefore we fit the distributions (examples are plotted as 
solid blue histograms in figure 11 with equation © using 
/ri, /r 2 , (Ti, (T 2 and ai /2 as free parameters. We do this for all 
mass-redshift bins given in table[T] At this point we are seek¬ 
ing the best fit for each bin individually. Examples of the re¬ 
sulting fits are overplotted in figureHas solid red areas. The 
values of the parameters for all of those fits are presented 
in figure [5] as a function of halo mass Mvir and redshift z. 
An analogue plot for the baryons was also produced, since it 


looks very similar to the gas plot it was omitted. The plots 
show that /ri as well as ^2 increase with increasing redshift, 
they also increase with increasing halo mass. A physical in¬ 
terpretation of these trends is that there is more inflow into 
heavier haloes and that there is also more inflow at earlier 
cosmic times. The Arp values for fii and fj .2 are considerably 
higher than their ART or MN counterparts. This is an effect 
of recycling. For the ci or (T 2 values the data are noisy and 
no clear trends are visible. Some of the values for ai /2 are 
also noisy. For gas however one can say that in more than 
three quarter of all the cases (17/22) ai /2 is above 0.75. The 
mean value for ai /2 averaged over all 22 mass-redshift bins 
from the three different suites is < ai /2 >= 0.82. These two 
numbers indicate according to the physical interpretation of 
the double-Gaussian curve that at least three quarter (pos¬ 
sibly rather four fifth) of the inflowing material is coming in 
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via smooth streams and only less than one quarter (rather 
one fifth) is coming in via merger events. This is^n agree¬ 
ment with the findings of iRomano-Di'az et alJ (|2014h who 
demonstrated that galaxies grow predominantly by smooth 
accretion from cosmological filaments. 

From the values shown in figure [5] we deduce that the 
mean accretion rate for the mergers (given by ^2) is larger 
than that of the smooth component (given by /ii). In log- 
space it is higher by ~ 0.5 which means that it is roughly 
a factor 2 — 3 larger. These values indicate in the terms of 
galaxy formation that significantly more material is coming 
into the host halo via smooth accretion than via merger 
events. The standard deviation of the merger accretion rate 
((J2) is 0.2 —0.3 dex, showing no trend with mass or redshift. 
For the smooth component (ai) it is 0.12 — 0.24 dex. There 
are only marginal differences between gas and baryons. The 
gas seems to account for the majority of the mass of the 
baryons. We calculated the gas fraction from the simulations 
and find that /gas = 85 %. This is the case for all three suites 
of simulations, MN, ART as well as Arp. 

We want to have a more general model and found that 
the behaviour with mass and redshift of two out of the five 
parameters, namely of fii and of p2 can best be described 
by a logarithmic relation: 

X = Ax ln (2 + l) + Bx ln(^^^ -bCx (4) 

It is valid for the two parameters and fi 2 , denoted by X, 
whereas Ax, Bx and Cx are the three free parameter de¬ 
scribing the logarithmic relation for both parameters. When 
fitting this equation we had to account for the very strong 
signs of recycling found in the Arp suite of simulations, as 
we already did for figure [3l There we showed the amount 
of inflow and found that the values for the Arp simulations 
had to be scaled down by some factor due to the effects of 
rescaling. In the following we are going to fit two of the pa¬ 
rameters of equation © which are also prone to be affected 
by recycling. The positions of the centres of the first and 
the second Gaussian of the distribution on the inflow (pi 
and p2 as defined by equation [3]) in the Arp simulations are 
considerably higher than in the MN or the ART simulations. 
However the reader should note that those are completely 
different physical values than the average amount of inflow 
< M > which are related but their relation is non-linear and 
therefore the effect of recycling has a differently strong effect 
on each of them. I.e. we have to scale the values for < M >, 
pi and p2 for the Arp simulations by independent scaling 
factors. The pi values of the Arp simulation are a factor 
of /rec(pi) = 1.45 higher than the corresponding values for 
pi of the MN or the ART simulations. The p2 values of the 
Arp simulation are a factor of /rec(p2) = 1.20 higher than 
the corresponding values for p2 of the MN or the ART sim¬ 
ulations. Again both factors remain constant over the whole 
redshift and host halo mass range considered. We present 
the best-fitting values for Ax, Bx, Cx and /rec(pi/2) in 
table m 

In figure we compare the general model of equation 
© with the parameters of the double-Gaussians of the sim¬ 
ulated distributions. Show are the positions of the centres of 
the first (pi, upper panels) and the second (p2, lower pan¬ 
els) Gaussian of the distributions of the inflow as a function 
of mass (left panels) and as a function of redshift (right pan- 


X 

Ax 

Bx 

Cx 

/rec(/J'l/2) 

Mlgas 

0.87 

0.44 

-11.4 

1.45 

/^2gas 

0.79 

0.39 

-9.4 

1.20 


Table 2. The best-fitting values for two (pi and P 2 ) out of the 
five parameters of equation 13 when fitted with our model (equa¬ 
tion 11 to the data of all bins as given in tabled as shown in in 
examples in figure [4] (red solid areas and labels). The other three 
parameters (ui, cr 2 and ai/ 2 ) were far too noisy to allow any 
meaningful fits. The column /frec(/^i/ 2 ) indicates the factor the 
Arp suite of simulations was scaled down by to account for its 
strong recycling. The actual comparison of the model (equation 
0 using the tabulated values and the simulations can be seen in 
figure [fi] the accuracy is exceptional. 

els). There are very strong logarithmic dependencies present 
in all panels as described by our general model which is over- 
lain as solid black lines. The data points within each panel 
were scaled with the help of the scaling relation in the corre¬ 
sponding other panel to the values indicated. The colour bar 
axes indicate the values the data points used to have before 
rescaling. Please refer to figure [3 for the unrescaled data. 
The panels show that our simple phenomenological model 
(equation |4] and table [2j is an excellent description of the 
data. 


5 CONCLUSIONS 

In this paper we looked at the inflow rates of accretion along 
streams from the cosmic web into galaxies at high redshifts 
using three sets of AMR hydro-cosmological simulations. We 
calculated the amount of inflow as a function of radius, host 
halo mass and redshift. We computed the distribution of 
inflow rates and we found the following: 

• The inflow rates are roughly constant with radius, their 
behaviour with h ost halo mass and redsh ift follows clearly 
the predictions of iNeistein fc Dekell (l2008l l. 

• The distributions of the log accretion rates can be 
very well described by a “double-Gaussian” functional form 
(equation [3 that is the sum of two Gaussians, the pri¬ 
mary corresponding to “smooth” inflow and the secondary 
to “mergers”. 

• Most of the amount of inflow (> 80%) is entering the 
halo at low inflow rates (i.e. smooth accretion) and only a 
small portion (< 20%) of the inflow is coming in at high 
inflow rates (massive merger events). 

• Two out of the five parameters of the double-Gaussian 
function (namely /ri and ^ 2 ) have a strong redshift as well 
as host halo mass dependence. In equation (|4]), tableland 
figure |6] we present a simple phenomenological model, which 
describes the shape of the double-Gaussian distribution as 
a function of mass and redshift. 

• The standard deviation of the total accretion rate is 
0.2 — 0.3 dex, showing no trend with mass or redshift. 

• The suite of simulations that include strong feedback 
shows massive signs of accretion due to recycling: i.e. for 
the Arp simulations we see (compared to the MN or ART 
simulations that have only weak feedback) an increase in 
the amount of inflow by a factor ~ 4.0, independent of mass 
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and redshift. Due to the same effect an increase in fii of the 
order of ~ 1.45 and an increase in fi 2 of the order of ~ 1.2 
(see table EJ is also seen. 

• There are hardly any differences between the gas and 
the baryons. Gas seems to account for the majority of the 
mass of the baryons anyway and additionally the star for¬ 
mation seem to locally follow the gas. 

• There is more absolute inflow into heavier haloes com¬ 
pared to lighter haloes also there is more absolute inflow at 
earlier cosmic times compared to later cosmic times. 

A double-Gaussian distribution has already been found 
bv ISargent et ah] (l2012| j. who looked at the distributions of 
specific star formation rates at fixed stellar mass in obser¬ 
vations. They interpreted the bimodality of this functional 
form as contributions from either main-sequence star forma¬ 
tion or star formation during starburst activity. The analogy 
between their and our findings indicates that connections 
might exist between smooth accretion and main-sequence 
star formation on the one hand as well as accretion through 
merger events and starburst activity star formation on the 
other hand. However the reader should note that just having 
the same functional form is no proof. 

In our most advanced simulations, the Arp simulations, 
which include the strongest feedback recipes, we see strong 
evidence for the reaccretion of formerly ejected material. 
This suite of simulations has more massive outflows and part 
of the outflowing material rains back again to the central 
galax y. An effect coined “recycling” l|ODDenheimer et al.l 
l20ldj . The inflow in the Arp suite of simulations is enhanced 
due to recycling compared to the MN or the ART simula¬ 
tions. This effect provides a third distinct acc retion mode 
along with the cold and hot modes d escribed bv iKeres et al.l 
(|200^ and iDekel fc BirnboimI ll200dj . 

One potential limitation of our simulations may arise 
from the artificial pressure floor imposed in order to prop¬ 
erly resolve the Jeans mass. This may have an effect on the 
temperature and density of the densest and coldest parts of 
the streams, with potential implications on the estimated 
inflow rates. Still modern hydrodynamical simulation codes 
are the best available tools for recovering the stream prop¬ 
erties. With 17-70 pc resolution, and with proper cooling 
below 10'*K, these simulations provide the most reliable de¬ 
scription of the cold streams so far. 

In a fort hcoming companion paper (iGoerdt et al.l 
lin preparatiorJl we are going to address the role of merg¬ 
ers versus smooth flows by analysing the dumpiness of the 
gas streams. We evaluate each clump mass and estimate a 
mass ratio for the expected merger. Finally we will look at 
the distribution of constituents amongst the clumps. 

We conclude that gas is flowing into a galaxy’s halo 
mainly as smooth accretion flows with only a minority 
(< 25%) coming in via merger events such as clumps or 
small satellite galaxies. The accretion rate distributions can 
be described by a double-Gaussian decomposition (equation 
0. That this functional form is also used to describe the spe¬ 
cific star formation rate distributions in observations could 
indicate that it is the smooth accretion that fuels main- 
sequence star formation on the one hand and the accretion 
of material through merger events on the other hand that 
fuels star formation of starburst events. 
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